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Abstract 

Constraints on the Equation of State for symmetric matter (equal neutron and proton numbers) 
have been extracted from energetic colhsions of heavy ions over a range of energies. Cohisions 
of neutron-deficient and neutron-rich heavy ions now provide initial constraints on the EoS of 
neutron-rich matter at sub-saturation densities from isospin diffusions and neutron proton ratios. 
This article reviews the experimental constraints on the density dependence of Symmetry Energy 
at sub-saturation density. 

1 Introduction 



> 

^ The Equation of State (EoS) of cold nuclear matter can be written as the sum of the energy per nucleon 
of symmetric matter and a symmetry energy term, Es 

E{p,6) = Eo{p,6 = 0) + Es (1) 

where 6 = {pn — Pp)/p is the asymmetry; p„, pp and p are the neutron, proton and nucleon densities, 
respectively. The first term on the RHS, Eq{p,6 = 0), is the EoS term for symmetric nuclear matter 
with equal fractions of neutrons and protons. Significant constraints on the symmetric matter EoS at 
1 < p/po ^ 4.5 have been obtained from measurements of collective flow [1] and Kaon production [2]. 
However, well-determined constraints on the symmetry term, Eg = S{p, 5^), are few. S{p) describes the 
density dependence of Es- 

Since large variations in nuclear density can be attained momentarily in nuclear collisions, constraints 
on the Equation of State (EoS) can be obtained by comparing measurements to transport calculations 
of such collisions. The symmetry energy has been recently probed at sub-saturation densities via 
isospin diffusion [31 H], and by double ratios involving neutron and proton energy spectra [5j. These two 
observables largely reflect the transport of nucleons under the combined influence of the mean fields and 



1 



the collisions induced by residual interactions; thus, they should be within the predictive capabilities of 
transport theory. 

This article focuses on investigations of the symmetry energy at sub-saturation density. The con- 
straints obtained in heavy ion collisions are compared to constraints obtained from nuclear structure 
studies, from Isobaric Analog States, Sn isotope skins. Pygmy Dipole Resonance experiments, and the 
Giant Dipole Resonance data. I will briefly discuss investigations of the density dependence of symmetry 
energy at supranormal density. 



2 ImQMD: Transport calculations for heavy ion collisions 

Transport models have been used to describe the dynamics of heavy ion collisions. We chose to use the 
ImQMD code mainly because of its ability to describe fragments which are quite important in describing 
isospin observables. In the QMD model, nucleons are represented by Gaussian wavepackets. The mean 
fields acting on these wavepackets are derived from an energy functional with the potential energy U 
that includes the full Skyrme potential energy with only the spin-orbit term omitted: 

U = Up + Umd + Ucouh Up.md = j Up^mdd^^ (2) 

where, Ucoui is the Coulomb energy and Up^rad represents the nuclear contributions in local form and 

^ P . P P^^ , 9sur /j-j ^2 I 9sur,iso /-r-, /■ \\2 , , P n^, r2 , P ^ /o\ 

'''^2^^^^^2^^^P^ +— — (V(p„-Pp)) +Y^-r^ P + 9pr^ (3) 
Po V + ^ Po -^Po Po Po po' 

where the asymmetry 6, and p„ and Pp are the neutron and proton densities, respectively. A symmetry 
kinetic energy density of the form %^(^)^''^5^p and symmetry potential energy density of the form 



^^{—y6^p were used in the following transport model comparisons. The energy density associated 



with the mean-field momentum dependence is represented by 
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where are nucleon Wigner functions, Ap = \pi — P2\, the energy is in MeV and momenta are in 
MeV/c. The resulting interaction between wavepackets is described in Ref. [6]. Unless otherwise noted, 
we use a=-356 MeV, /3=303 MeV and r] =7/6, corresponding to a isoscalar compressibility constant of 
K=200 MeV, and ^,^^=19.47 MeVfm^, ^,™=-11.35 MeVfm^, C,,fc=24.9 MeV, C,^p=35.19 MeV, and 
gpr =0 MeV. 

These calculations use isospin-dependent in-medium nucleon-nucleon scattering cross sections in 
the collision term and Pauli blocking effects that are described in [71 [8]. Cluster yields are calculated 
by means of the coalescence model widely used in QMD calculations in which particles with relative 
momenta smaller than Pq and relative distances smaller than Rq are coalesced into one cluster. In the 
present work, values of Rq = 3.5 fm and Pq = 250MeV/c are employed. 



3 Experimental constraints from heavy ion collisions 

We turn our attention first to the interpretation of neutron/proton double ratio data, which derives 
its sensitivity to the symmetry energy from the opposite sign of the symmetry force for neutrons as 
compared to protons [9] . First experimental comparisons of neutron and proton spectra in Ref. [5] used 
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a double ratio in order to reduce sensitivity to uncertainties in the neutron detection efficiencies and to 
relative uncertainties in energy calibrations of neutrons and protons. This double ratio, 



DR{nlp) = R^/p{A)/R^/p{B) 



dMn{A)/dE, 
dMp{A) /dE, 



'cm. 



'cm. 



dMn{B)/dE, 
dMp{B)/dE, 



'cm. 



'c .m. 



(5) 



was constructed by measuring the energy spectra, dM/dEc.m., of neutrons and protons for two systems 
A and B characterized by different isospin asymmetries. The star symbols in the left panel of Figure 1 
show the double ratios measured at 70° < 9cm < 110° as a function of nucleon center-of-mass (cm.) 
energy for central collisions of ^'^^Sn +^^^ Sn and ^^"^Sn +^^^ Sn at E/A=50 MeV. 

We have performed calculations at an impact parameter of b=l, 2 and 3 fm at an incident energy of 
50 MeV per nucleon for two systems: A=^^'^S'n +^^^ Sn and B=^^'^Sn +^^^ Sn. About 60,000 events are 
simulated for each impact parameter. Within the statistical uncertainties, the double ratio observable, 
DR(n/p), is nearly independent of the impact parameters within the range of 0<b<5 fm. The lines in 
the left panel of Figure 1 show the predicted double ratios DR(np)=i?„/p(124)/i?n/p(112) as a function 
of the cm. energy of nucleons emitted at 70° < 9c. m. < 110° for 7^=0.35, 0.5, 0.75, 1 and 2. The 
uncertainties for these calculations in Figure 1 are statistical. As fewer nucleons are emitted at high 
energy, the uncertainties increase with increasing energy. Despite the large experimental uncertainties, 
especially for the higher energy (>40 MeV) data, the trends and magnitudes of the data points definitely 
rule out the very soft ( 7^=0.35, dotted line with open diamond points) and very stiff ( 7^=2, dotted line 
with closed diamond symbols) density-dependent symmetry terms. We find the 2a uncertainty range 
of values for 7^ to be 0.4< 7^ < 1.04, corresponding to an increase in by 4 above its minimum of 

~ 2.1 near 7^ = 0.7. 

The density dependence of the symmetry energy can also be probed in peripheral collisions between 
two nuclei with different isospin asymmetries by examining the diffusion of neutrons and protons across 
the neck that joins the nuclei. This "isospin diffusion" generally continues until the two nuclei separate or 
until the chemical potentials for neutrons and protons in both nuclei become equal. To isolate diffusion 
effects from other effects such as pre-equilibrium emission. Coulomb effects, and secondary decays, 
measurements of isospin diffusion compare "mixed" collisions involving a neutron-rich nucleus A and 
a neutron-deficient nucleus B to the "symmetric" collisions involving A+A and B+B. The degree of 
isospin equilibration in such collisions can be quantified by rescaling the isospin observable X according 
to the isospin transport ratio Ri{X) [3| given by 



In the absence of isospin diffusion, Ri{XA-\-B) = Ri{XA+A)=^ for a collision involving a neutron-rich 
projectile on a proton-rich target. Likewise, without diffusion, Ri{XB+B) = -Ri(-^B+A)=-1- On the 
other hand, if isospin equilibrium is achieved for roughly equal sized projectiles and target nuclei, 
Ri{XA+B) = Ri{XB+A) ~ 0. Using ImQMD model, we investigate Ri{6) as a function of the rapidity 
of the emitted particles. The results of the calculations are shown by the lines in the right panel of 
Figure 2. The data, represented by the star symbols, are obtained from the yield ratios of the mirror 
nuclei pair, '^Li and ^Be, Rj = Ri{X = ln(Y(^Li)/Y(^Be)))^. This calculation reproduces the shapes 
and magnitudes of the rapidity dependence of the isospin transport ratios Rj. The analysis brackets 
7j values to lie in the region 0.4< 7^ <1.0 using the same 2a criterion. There is substantial overlap 
between the constraints obtained from isospin diffusion and double neutron and proton yield ratios. 

Constraints on the exponent ■ji depend on the symmetry energy at saturation density, Sq = S{po) = 
0.5 * {Cs^k + Cs^k)- Increasing Sq has the same effect on the isospin transport ratio as decreasing 7j. 
Thus 5*0 and 7^ cannot be determined independently. To provide constraints on 5*0 and 7^ parameters, 
contours with constant values can be created by doing two-dimensional analysis in the 7^ — 5*0 



Ri{X) = 2 



X — {Xa+a + Xb+b)/'^ 
Xa+a — Xb+b 



(6) 
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Figure 1: (Color online): Left panel: Comparison of experimental double neutron-proton 
ratios (star symbols), as a function of nucleon center of mass energy, to ImQMD calculations 
(lines) with different density dependence of the symmetry energy parameterized in terms of 
7j in Eq. (1). Right panel: Comparison of experimental isospin transport ratios obtained 
from the yield of A=7 isotopes (star symbols), as a funtion fo the beam rapidity, to ratios 
from ImQMD calculations (lines) for b=6 fm. 



parameter space. Unfortunately, the CPU requirements for such large-scale ImQMD calculations are too 
intensive to make such analysis practical. Instead, We have preformed a series of ImQMD calculations 
at b=6 fm with different values of 7^ over selected values of Sq between 25 to 40 MeV to locate the 
approximate boundaries in the 5*0 and 7^ plane that satisfy the 2a criterion in the analysis of the 
isospin diffusion data. The shaded area shown in the inset of Figure 3 is the result of such analysis. 



4 Comparison of experimental constraints at sub-saturation 
density 

Near saturation density, one may expand the symmetry energy, S{p), about the saturation density, po, 

S[p) = Sq + - + —r^y +••• 7 

where L and Kgym are slope and curvature parameters at po. The slope parameter, L, is related to po, 
the pressure from the symmetry energy for pure neutron matter at saturation density: 

L = 3po\dS{p)/dp\po = [3/po] ■ po (8) 

Po, provides the baryonic contribution to the pressure in neutron stars [lO], where the energy of sym- 
metric matter, Eq{p,S = 0), does not contribute to the pressure. It is also related to the neutron skin 
thickness ( SRnp) of neutron-rich heavy nuclei including ^osp^ [T2] . 
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Figure 2: (Color online): Representation of the constraints on the density dependence of 
symmetry energy as a function of the symmetry energy at saturation density in neutron 
matter. The hatched areas are extracted from recent analysis of nuclear collisions of Sn 
nuclei. The two dashed lines represent the constraint boundaries obtained from the analysis 
of the skin thickness of tin isotopes [13]. The horizontal line indicates the nearly constant 
mean L value in the overlap region. 



The series of ImQMD results calculated at b=6 fm with different values of 7^ over selected values of 
5*0 between 25 to 40 MeV give the approximate boundaries in the So and L plane that satisfy the 2a 
criterion in the analysis of the isospin diffusion data. The area bounded by the two solid diagonal 
lines in Figure 2 and the shaded region in Figure 3 represents a conservative estimate in such effort. 
Another independent constraint will be needed to provide constraints in Sq or L. 

Recent analysis of the skin thickness of tin isotopes [I3] suggests that L decreases with increasing 
5*0. The two diagonal dashed lines plotted in Figure 2 represent a 2cr analysis of the discussed in 
Ref. p^. Since the trend is opposite to that of the constraints provided by the heavy ion reactions 
(solid lines), the overlap of the two constraints provide a range of Sq and L values. As pointed out by 
Chen et al, the interesting thing to note is that the mean L values, represented by the horizontal line 
in Figure 2 obtained from the overlap region of the two constraints, are nearly constant. The range of 
L and 5*0 values in the overlap region corresponds to a mean skin thickness of '^^^Ph of about 0.18±0.04 
fm [TU]. The range of experimental neutron skin thicknesses of "^^^Ph varies from about 0.13 to 0.22 fm 
[m [151 [m 113 HH] which are consistent with results predicted by the constraints. 

We have included, in Figure 3, other constraints from the properties of individual nuclei. The 
corresponding values of pure neutron matter pressure, po? are given in the right axis. The lower box 
centered at 5*0 = 32.5 MeV depicts the range of Sq and po values from analysis of low-lying electric 
dipole strength (Pygmy Dipole Resonance, PDR) in neutron-rich nuclei ^^Ni and ^^"^Sn [121 ISO]- The 
upper box centered at Sq = 32.5 MeV depicts the constraints reported in Ref. [21] from the analysis of 
the excitation energies of Isobaric Analog states (IAS). All these constraints overlap significantly with 
the hatched regions obtained from the heavy ion collisions of tin isotopes as described in section 3. The 
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Figure 3: (Color online): Comparison of various constraints on the density dependence of 
symmetry energy (inset) and pressure as a function of the symmetry energy at saturation 
density in neutron matter. The hatched areas are extracted from recent analyses of nuclear 
collisions of Sn nuclei. Open rectangles and the symbol in the inset were obtained from the 
properties of individual nuclei. The inset shows the density dependence of the symmetry 
energy of the shaded region. The symbol in the inset represents the GDR results |22j . 



inset shows what this analysis implies about the density dependence of the symmetry energy. A value 
for the symmetry energy at around 0.6 normal nuclear matter density, extracted from Giant Dipole 
Resonance (GDR) data [22], is indicated by the symbol in the inset. 



5 Summary and outlook 

In summary, the ability of QMD models to reproduce both isospin diffusion data and double ratio of 
neutron and proton spectra data is an important step forward in obtaining information about the sym- 
metry energy in heavy-ion collisions. The heavy ion as well as nuclear structure observables examined 
here provide consistent constraints on the density dependence of the symmetry energy. 

To explore the density dependence of the symmetry energy, an international collaboration, consisting 
of an interdisciplinary, experimental and theoretical team of scientists, has been formed [23j to conduct 
a series of experiments at unique facilities based in the United States (the National Superconducting 
Cyclotron Laboratory at Michigan State University), Japan (the Radioactive Ion Beam Factory at 
RIKEN) and GSI, Germany. Each facility enables the exploration of a different density range. With 
experiments underway at GSI, better constraints at supra normal density is expected in the near future. 
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